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Abstract 

We analyze the operation of the wireless single-electron logic family based 
on single-electron-parametron cells. Parameter margins, energy dissipation, 
and the error probability are calculated using the orthodox theory of single- 
electron tunneling. Circuits of this family enable quasi-reversible computation 
with energy dissipation per bit much lower than the thermal energy, and 
hence may circumvent one of the main obstacles faced by ultradense three- 
dimensional integrated digital circuits. 
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I. INTRODUCTION 



Single-electron tunnelingm in systems of small-capacitance tunnel junctions has at- 
tracted much attention during the last decade because of both the interesting physics and 
the possibility of useful devices. Today the basic physics of single-electron tunneling is suffi- 
ciently well understood, and its possible applications (for reviews see, e.g., Refs.i and@) are 
becoming an important issue. The practical value of several analog applications has already 
been proven: for example, the single-electron transistori'i as a highly sensitive electrometer ,0 
the single-electron pumpi as a standard of dc current,! and the array of tunnel junctionjisl 
as an absolute thermometer.0 

In the field of digital electronics, however, only rudimentary devices have been demon- 
strated so far.llHli The recent inventionEl of hybrid circuits based on single-electron- 
transistor readout from fioating-gate memory cells presumably opens a way to room- 
temperature memories with density beyond 10^^ bit/cm^. (This concept has been supported 
by recent experimentslll with the first low-temperature prototypes of such memory cells). 
The potential of single-electronics in logic circuits, however, remains uncertain due to sub- 
stantial problems which have to be solved for both types of single-electron logic devices 
suggested so far. 

Devices of the first type use single-electron transistors instead of field effect transistors in 
"voltage-state" circuits similar to conventional digital electronics.0 In these devices the effects 
of charge quantization are restricted to the interior of the single-electron transistors. As a 
result, the design of such circuits (including complementary!^) is relatively straightforward. 
They suffer, however, from a relatively high static power consumption. 

The second possible way (which allows to reduce the power dissipation), is to code bits 
directly by single electrons everywhere in the circuit .iHH (A different, interesting approach 
based on rf parametric excitation of single-electron-tunneling oscillations has been suggested 
by Kiel and Oshima;!! unfortunately their original proposal seems to have run into consid- 
erable implementation problems^). It is important that these "charge-state" logic circuits 
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do not necessarily require long wires and can be implemented using only small conducting 
islands separated by tunnel barriers, with the necessary power and timing provided by an 
external rf electric field.0 Recently we suggested^ a new device for charge-state, wireless 
logic circuits - the Single Electron Tunneling Parametron (or SET Parametron). In com- 
parison with the wireless single-electron devices suggested earlier,^ the SET Parametron 
may have wider parameter margins and extremely low power dissipation, lower than the 
thermodynamic "limit" of ksT \n2 per bitS 

In this paper we present the results of a detailed analysis of the SET Parametron and 
basic logic circuits using this device. The major characteristics of the SET Parametron, 
including its parameter margins, speed, error rate, and power consumption, will be discussed. 



II. SET PARAMETRON: THE BASIC IDEA 

The basic cell of the new logic family, the SET Parametron, consists of three conducting 
islands - see Fig. |l|a. (In our numerical calculations we have assumed the islands to be of 
spherical shape, with R being the radius of the edge islands 1 and 3, and r the radius of the 
middle island 2; this assumption is, however, not important for the device operation.) The 
middle island is slightly shifted off the line passing through the centers of the edge islands. 
We will identify this line with x-axis, and the direction of the middle island shift with y-axis 
(which is the symmetry axis of the cell). Electrons can tunnel through small gaps between 
the middle and edge islands, but not directly between the edge islands because of their much 
larger spacing (Fig. |l|a). 

Let us assume that the cell as a whole is charged by one uncompensated electron. (This 
assumption makes the explanation of the operation principle simpler; later we will consider 
the more natural case of an initially neutral cell.) If the cell is biased with a sufficiently 
strong "clock" electric field Ec > directed along axis y, the electron is obviously located at 
the central island.S Now let the field be decreased gradually so that eventually it changes 
direction to negative. At some moment the electron will have to tunnel to one of the edge 
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islands because these states become energetically preferable. Because of the geometrical 
symmetry, the choice of island (left or right) will be random, i.e. the charge symmetry of 
the system will be broken spontaneously. 

However, if there is a weak "signal" field Eg along direction x (this field may be applied, 
for example, by the neighboring similar cell), the final position of the electron will depend on 
the sign of Eg. A natural way to discuss this effect is to say that the signal field Eg creates 
an energy difference A between the electron states in islands 1 and 3, and that the electron 
prefers to tunnel into the island with the lowest energy state (Fig. |l|b). If now the clock field 
Ec becomes negative and large, it creates a high energy barrier W (see Fig. [l|b) between 
the edge islands, so that the electron is reliably trapped in the island it has tunneled to, 
regardless of further changes of the signal field Eg. 

Thus if fluctuations in the system are low enough, and the clock field changes slowly 
enough, even a small field Eg of the proper sign at the decision-making moment (when 
W{t)=0) is sufficient to ensure a certain robust final polarization of the cell. This process 
can be interpreted as a reliable recording of one bit of information; for example, the electron 
on the right island may mean logical "1" while the electron on the left island encode logical 
"0". Now the dipole moment of the cell can in turn be used to produce the signal field 
to control the other cells during their decision-making moments, and hence determine their 
information contents (see Section V below). 

This "parametric" amplification of signal, which gave the SET Parametron its name, is 
quite similar to the operation of a Josephson junction device called Parametric Quantronim 
(re-invented later as the "Flux Quantum Parametron" S), and in a broader sense to rf-driven 
parametrons.iiS The main difference is that all the previously analyzed parametrons are 
described by a continuous degree of freedom (e.g., the Josephson phase (p in the Parametric 
Quantron) . This variable can take any value, and the discrete information states correspond 
to energy minima on the (f) axis. On the other hand, in the SET Parametron with low tunnel 
conductances, the possible states are discrete (the electron has to be definitely on one of the 
islands at each particular instant). 



III. SET PARAMETRON: QUANTITATIVE ANALYSIS 

Let us consider the operation of one SET Parametron cell using the orthodox theory0 
which is adequate for the description of single-electron tunneling in systems with many 
background electrons in each island and sufficiently low conductance G of tunnel junctions 
{G <^ e'^/h). In this theory the electron is always localized in one of the islands, and the rate 
r of each tunneling event is solely determined by the corresponding decrease W = —Ae of 
electrostatic energy e of the system. Using elementary electrostatics, we obtain the following 
expressions for the energy decreases corresponding to four possible tunneling events in the 
SET Parametron (tunneling to/from each of the edge islands is denoted with signs + and 
— , respectively): 

W^± = ±e(02-03)-lV23, (1) 

4>i = j2iC-%q,-ED,, (2) 



e2 

Wii = — 
u 2 



{C-% + {C-'),, - 2{C-\] , (3) 



Qi = Qio + . (4) 

Here l{r) means left (right) junction, 0j is the electrostatic potential of i-th island (for 
numbering, see Fig. |l]a), W12 and W23 are the Coulomb blockade energies, and is the 
inverse capacitance matrix. Vectors eDi = {eDix, eDiy, eDiz) are the dipole moments of 
islands when charged with single electrons, E is the external electric field (accepted to be 
uniform) and Qi is the total electric charge of i-th island, while g^o is the uncompensated 
part of its background charge. 

In the limit when the island radii are smaller than the spacing between islands, the 
following simple approximation for the capacitance matrix elements can be used: 
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A = r,, (5) 

where = (x,, Zi) is the i-th island center, and e is the dielectric constant of the dielectric 
environment. In the case of comparable radii and island spacing, Cij have to be calculated 
numerically, for example, using the method of multiple electrostatic images (see Appendix 
A). The results show that approximation gives accuracy better than 10% if the spacing 
between the islands is larger than the largest island's radius. 

In our initial analysis (until Section VII) we will assume that temperature T is low and 
that the electric fields are changing slowly (adiabatically). In this case the full result of the 
orthodox theory, 

^ GW 

e2 [1 - exp(-iy/A;BT)] ^ 

is reduced to the simple rule that the system always follows the local minimum of the total 
electrostatic energy. 

Figure |^a shows the phase diagram of stationary charge states of a cell charged by 
one uncompensated electron. Each state is locally stable within a diamond-shape region 
corresponding to two conditions: 

|e(02-0i)| < W^12, |e(02-03)| < W^23, (7) 

being valid simultaneously, so that no tunneling event is possible, W^^ < 0. Within the 
small-island approximation (^ the position of the diamond center is 

111/ N / 1 1 

Es = -Aril — nsje — 

' 47reeo 2 ^ ^' \2d R 

1 1 1 



47reeo 2 h 

/I 1 

+2^26 

\r s 



N /2 1 1 



(8) 



while the diamond width and height are 

e a e 

(here s = Si2 = S23 = {cP + 6^)^/^ is the distance between island centers - see Fig. ^a). 
The diamonds are periodic along both axes Eg and Ec-, with periods 

11/1 1 \ 

= : e , 
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Equations (||)-(|10D are valid even if the approximation (|) is not used; however, in this 
case l/R, 1/r, 1/s, l/2d (divided by 47reeo) need to be replaced by Cii, C22, C12: 
respectively, and also b and d by D2y — Diy and D2X — Di^ (the corresponding equations are 
direct consequences of Eqs. (P-(0) and (|^)). 

The period along axis Eg corresponds to the transformation (rii, n2, n^) — > (rii — 1, n2, ^3+ 
1), while the period along axis corresponds to (ni,n2,n3) (711 — 1,^2 + 2,^3 — 1). Notice 
that these transformations allow only even-number changes at the middle island; the "odd" 
set of diamonds can be obtained from the "even" set by the shift {AEs/2, AEc/2). 

Figure l^a shows the results of more exact, numerical calculations of the phase diagram for 
a particular cell; however, they are quite close to the approximate analytical results (^-(|T0|). 
(The approximation accuracy is better than 10% and it becomes even better than 1% after 
the normalization by AEg and AE^..) For example, the minor offset of the diamond corners 
from the horizontal axis, which arises due to deviations from the analytical approximation, 
is hardly visible. 

The most important qualitative feature of the phase diagram is the existence of bistability 
regions (shaded in Fig. ^) where either of two charge states is locally stable. We will refer 
to the SET Parametron in these regions as being in an ON state, while the remaining part 
of the phase diagram corresponds to one of the possible OFF states of the system. If signal 
field Eg is low (as we suppose in the following), only the set of diamonds along axis Ec 



may be implemented. The arrow in Fig. ^a shows a possible way (described qualitatively 
in Sec. II above) of switching the cell from a monostable OFF state to a bistable ON state 
by changing the clock field Ec- It is clear that the sign of a small Eg field determines which 
diamond boundary will be crossed first, and hence determines the charge state of the system 
within this bistable region. 

Notice that crossing any phase boundary outside of the bistable regions involves one 
tunneling event with negligible energy difference W. In the adiabatic limit this transition 
is reversible, and the associated energy loss is infinitesimal (see Sec. VII). However, leaving 
the bistable state from the state with higher energy results in two sequential tunneling 
events and leads to a finite energy loss. This irreversible switching should be avoided in 
reversible computation. Below we will show that SET Parametron circuits can operate 
using exclusively reversible transitions. 

IV. SINGLE EXCITON PARAMETRON 

While the case of the cell charged by one electron is the simplest for understanding the de- 
vice operation, a more natural option is the initially electrically neutral cell: 7^1+^2+^3 = 0. 
(This variety of SET Parametrons may be called the "Single Exciton Parametron", since 
digital bits in it are presented by electrostatically bound electron-hole pairs, although the 
physics of this bound state is rather different from that of the usual excitons in semiconduc- 
tors.) Figure 03 shows the phase diagram of such a cell with the same geometrical parameters 
as in Fig. 0a. The only difference between the diagrams is a change in the charge states 
labeling and a fixed shift along Ec axis, so that the diagram in Fig. is symmetric about 
both axes. (Actually, any change of the initial charge of the cell is equivalent to a fixed 
offset of the fields Ec and Es). 

Due to its field sign symmetry, the Single Exciton Parametron allows two natural modes 
of operation. The first is illustrated by the arrowed rectangle in the center of the diagram 
(instead of the thick line in Fig. ^a we have drawn the rectangle in Fig. |^ to emphasize that 
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Es may shift during the operation). The periodic clock field Ec{t) causes two switchings 
OFF^ON and two switchings ON— >OFF per period. Another possible mode of operation 
is shown by the smaller rectangle centered at Es=AEs/2. It also provides two pairs of 
ON/OFF switchings per period of the clock field Ec{t), with the bistable state implemented 
at small fields. In our analysis we will concentrate on the former operation mode because it 
does not require the additional dc bias of the signal field, though the latter mode may allow 
wider parameter margins. 

Figures |^ and §3 show the phase diagrams of SET Parametrons for a particular set 
of geometric parameters. Changing the parameters may result not only in quantitative but 
also qualitative changes of the diagram. In particular, regions with three stable charge states 
appear if the ratio 6Es/AEs of the diamond width and period becomes larger than 2; if the 
ratio is larger than 3, four states can coexist, and so on. Geometrically the multi-stability 
requires a smaller middle island to increase the Coulomb barrier for the electron transfer 
between the edge islands via the middle island. For the sake of simplicity we will limit 
ourselves to the bistable case.0 



V. SHIFT REGISTER 

Figure |^ shows the standard scheme of operation of a shift register using parametron-type 
binary cells.i^ In this figure, each symbol represents a parametron cell at some particular 
phase of the periodic switching process. Dash means the monostable (OFF) state, while li 
means the i-th data bit which can be either "0" or "1" . Evolution of the symbols shows the 
bit fiow when the cells are periodically switched ON and OFF using periodic clock signals 
which are shifted by 27r/3 for each next stage of the register. 

In the initial state (upper line) each bit is encoded by the corresponding ON state of only 
one parametron. The signal field created by this cell is exerted on both neighboring cells, 
initially in their OFF state. Due to this field, when the right neighboring cell is switched 
from OFF to ON by the clock field (second line in Fig. |^), its logic state will be determined 
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by that of the initial celL Now when the information has been copied to the right cell, 
the initial cell can be switched OFF (third line of Fig. |^). Thus, the information has been 
moved by one cell (one register stage), and the process may be repeated again (lines 4 to 7) 
and again. Notice that unidirectionality of the bit propagation is achieved by the "running 
wave" of the clock while the structure itself is ID-isotropic. 

During each period of the three-phase clock, the information is shifted by three stages, 
so that we need 3 cells per bit of information. It is straightforward to increase the number 
of phases beyond 3. Typically this would make the system more robust, but require propor- 
tionally more hardware per bit. (More generally, the number of cells per bit has not to be 
integer: by changing the clock phase shift it can have any value larger than 2). 

Figure ^ shows a natural implementation of shift register using Single Exciton 
Parametron cells.ii Each next cell, the direction of the middle island shift is turned by 
7r/3 (within plane y — z). The system is driven by the clock field vector E{t) rotating in 
the same plane. This rotation provides the shift of E^. (which is now the component of the 
clock field in the plane of the corresponding cell) by 1/6 of the clock period for each next 
cell. When a cell is in ON state, its dipole electric moment creates the signal field Eg which 
is especially strong for its nearest neighbors, and thus determines the direction of electron 
tunneling when one of the neighbors in turned ON. 

The operation of the system is illustrated in Fig. ^. Three sine lines show the in-plane 
components of the clock field for each of three neighboring cells. Since the operation cycle 
corresponds to the central rectangle shown in Fig. 03, it differs slightly from the traditional 
mode shown in Fig. |^: each cell is switched OFF and ON twice per clock period.0 In the 
simplest case of weak coupling when the signal field Es is much smaller than the clock field 
amplitude, switchings happen when the clock field E^. crosses one of the thresholds -^5Ec/2 
and the cell enters or leaves the bistable region (shaded in Fig. For correct operation of 
the shift register, switching OFF— >0N of a cell (see, e.g., point 5* at the thicker middle curve) 
should occur when the previous neighboring cell is in ON state while the next neighboring 
cell is in OFF state, just as shown in Fig. |^. 
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In order to find the width of parameter window corresponding to correct operation of 
the parametron, let us first neglect the inter-cell interaction. Using the operation diagram 
shown in Fig. ^ and the phase diagram of Fig. |^ it is straightforward to obtain that the 
amplitude A of the rotating clock field has to satisfy the following three conditions: 

Asm( ) > 



2 2' 2 ' 
^sm(-)<— , 

A<Ai?.-^. (11) 

where 9 is the clock phase shift between the adjacent cells. The first condition is necessary 
to have the signal source cell still in ON state when a recipient cell switches from OFF to 
ON. The second relation ensures that at that instant the next cell is already in OFF state. 
Finally, the third inequality means that the clock field does not reach the next diamond at 
the Ec axis - see Fig. gb and Eqs. and (PD]), so that no undesirable charge states appear 
during the operation. 

Equations ( |lTD can be reduced to the following equation for the one-side margin ^ of the 

field amplitude, ^ = {Amax ~~ Amin) / i^Amax + Amin)'- 

, TT 9 2/\Ej5Ec - 1 " 1 / cos(^/e) 



^ = min 



tan(- 



(12) 



■4 2" 2AEe/5Ec - 1 + 1/ cos(^/e) 
Within the approximation d) the ratio I\EjbE^ is equal to 2 - (1/i? - 1 /2rf) /(1//2 + 1/r - 
2/s); this ratio increases when the radius of the middle island decreases. In the general 
case when the terms 1/-R, 1/r, l/2(i, and \ j s are replaced for the elements of the inverse 
capacitance matrix, the upper bound for /S.Ecl5Ec is equal to 2, while the lower bound is 
equal to 1.0 

For the particular case shown in Fig. ^ we have 9 = 27r/6, hence, Eq. (|I2[) gives the 
maximum margin of ^ = tan(7r/12) ^ 27% for the amplitude A of the rotating clock field if 
AEc/SEc > 1.5, so that the minimum in Eq. ( [T^ ) is determined by the first term. (In order 
to satisfy this additional condition the middle island should be sufficiently smaller than the 
edge islands.) 
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The above single-cell analysis is only an approximation, because in the real structure all 
the cells interact with each other and also the presence of neighboring cells somewhat changes 
the cell electrostatics. This is why we have carried out extensive numerical simulations of 
reasonably long shift registers (18 three-island cells). For a given geometry of the circuit 
we first calculated numerically the inverse capacitance matrix C^^^ and the vectors Di (see 
Appendix A); then these numbers were fed to a Monte-Carlo program which simulated single- 
electron tunneling events in the circuit using the "orthodox" theory.0 The logical input was 
a relatively long (typically, 16-bit) quasi-random bit sequence, repeated periodically many 
times. 

The simulations have shown that for the geometry presented in Fig. ^, elliptic polarization 
of the clock field (rather than the simple circular polarization implied above) is essential to 
provide broad margins. More exactly, the out-of-plane component should be larger than 
the in-plane component Ay, because of the field screening in z direction due to the finite 
size of the conducting islands. For example, for the circuit with r/R = 0.3, d/R = 1.5, 
b/R = 0.7, and the intercell period in y-direction l/R = 3, the best ratio A^/Ay is close 
to 1.5. If this is done, sufficiently wide margins for the amplitudes of the bias field are 
really achievable. For A^/Ay = 1.5, the 18-cell shift register with the parameters specified 
above operates correctly (at zero temperature and sufficiently low clock frequency) within 
the range 3.25 < Ay/ {e/AneeQR^) < 5.1, corresponding to a margin of ±22%. If now Ay is 
kept constant at Ay = 4.15 e/47reeo-R^, the allowed range for Az/{e/4:7ieeoR^) is between 4.6 
and 8.55, corresponding to a margin of ±30%. These margins are much broader than those 
for the first generation of wireless single-electron devices.0 

The margins for parameters which destroy the geometric symmetry of the parametron 
cells are much narrower. Numerical simulations show that for a shift register with the 
parameters quoted above, the maximum allowable shift of the edge islands of a cell in the 
opposite directions (parallel to the central island shift) is about ±0.010 i?. This corresponds 
to the maximal fiuctuation of ±0.4 degrees of the alignment of the line connecting the edge 
islands of the cell and the x-axis. (A shift of only one island leads to the same effect). 
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Possibly the situation may be improved by using alternative geometries and charge modes. 



VI. SET PARAMETRON LOGIC GATES 

The shift register described in the previous Section is actually a line of inverters. To 
have a complete set for the arbitrary logic functions we need to have other logic gates (e.g., 
NAND or NOR), and a circuit with a fan-out more than one ("splitter"). 

All these functions can be naturally implemented using the geometry shown in Fig. If 
the clock field rotation causes the signal propagation from the bottom to the top we get a 
fan-out-two circuit, because the dipole moment of cell F (in its ON state) will determine the 
charge state of both cell A and cell -B, set during their OFF— >0N switching. On the other 
hand, if the signal propagates from top to bottom, we get the implementation of a binary 
logic function (either F = A.NOR..B ot F = A.NAND.B, depending on the asymmetry 
provided by the background charges or dc bias field (Eg) imposed on cell F). 

In order to verify the operation of this system, we have carried out its numerical simula- 
tions in the symmetric single-exciton mode, using the same parameters as above {r/R= 0.3, 
d/R = 1.5, b/R = 0.7, l/R = 3) for each of three shift registers, each of them 6 cells long. 
The 2;-shift between upper registers was 6-R, while the y-spacing between cells A (B) and F 
was 2R, i.e. slightly less than the spacing I = 3R between cells in each register. 

The simulations have shown the following margins for the amplitude components Ay 
and Az of the clock field (at zero temperature and relatively low speed). For the fan-out- 
two gate, at Ez/Ey = 1.5, we get correct operation at 3.45 < Ay/^e/ATTeeoR"^) < 4.2, 
corresponding to a margin of ±9%. For the constant Ey = 4.15 e/47reeo-R^, the z-component 
should be within the range 4.6 < Az/{e/47TeeQR'^) < 6.3 (the margins of ±16%). For 
the NOR gate, with Ey = 4.1 e/47reeo-R^, E^ = 6.15 e/47reeo-R^ (this operating point is 
within the margins for both the shift register and fan-out circuits) correct operation is 
achieved for the dc bias field (Es) equivalent to the background charge shift in the range 
0.034 < qq/c < 0.053 (go = qo,ieft = —%,right)- For the NAND gate with similar parameters 
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the range is 0.035 < — go/e < 0.069. The difference between the results for the ranges 
(which should coincide because of the layout symmetry) is explained by the finite set of 
tested input signal sequences which did not have a symmetry between logical zero and unity. 
Hence, the true range cannot be wider than the common part of two ranges [0.035,0.053], 
and thus the real margin is slightly below ±20%. 

VII. ENERGY DISSIPATION DURING OFF^ON SWITCHING 

Let us calculate the energy dissipation during the switching OFF— i>ON of one SET 
Parametron cell in a fixed signal field Es- The switching is illustrated in Fig. ^ which shows 
the energies of three possible states of the cell as functions of time. For a cell charged 
by one electron the letters /, r, and m correspond to the location of the extra electron. 
For the single-exciton case we can use the same notation, counting one hole or electron as 
background charge. 

Within the framework of the "orthodox" theory the switching is described by the master 
equation,0 

d ^ 

— = cr^r^ + o-^r^ - crm(r/' + r+) , 

am + cri + ar = l, a„X-oo) = 1 , (13) 

where am, (^i, are the probabilities of finding the system in "m", and "r" states. The 
tunneling rates F^"^^^ are given by Eq. (|^), in which we now assume the linear dependence of 
the energies on time: 

W+{t) = -W-{t) = W{t) = at (14) 

W+it) = -Writ) = Wit)-A. (15) 

(This assumption is evidently valid only if the energy difference A due to the signal field 
in the moment of switching is much less than the maximal energy due to the rotating bias 
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field.) We also assume that the signal field makes state "r" energetically more preferable, 

= A > 0. 

If the clock speed and temperature are sufficiently low so that the system switches cor- 
rectly from "m" to "r" (bit errors will be considered in Section IX), then the energy dissi- 
pation can be calculated considering only states "r" and "m" and neglecting state "/". In 



this case the master equation (0) is simplified with the condition cri{t) = 0, 

d 



(1 - (T^)r+ - GrV- (16) 

(here we omit index r in tunneling rates for simplicity), and has an explicit solution 



J-oo . Jt' 



dt' , 



(17) 



where 



rs(t) = r+ + r- = r+ [1 + exp(-i^(t)/T)] 



The statistical average of energy loss by the moment t is given by the expression 



(18) 



W{t')&r{t')dt' 



(19) 



Using Eq. (|T7p we find, finally: 



l + exp(-^; 



l + exp(-m; 



xexp[-£ r^{t"')dt"']dt" 



dt'. 



(20) 



The solid line in Fig. ^ shows the total energy loss S = S{oo) as a function of the 
dimensionless switching speed 



(3 = ae^/GT^ 



(21) 



In the low-speed limit, /5 <^ 1, the switching process is almost adiabatic. It consists of 
numerous tunneling events (back and forth between "m" and "r" states) occurring within 
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the energy interval on the order of the thermal energy T around the point W{t) =0. (Figure 
1^ shows a particular Monte-Carlo realization for j3 = 0.1.) It is easy to check that in the 
purely adiabatic case when ar is determined by thermodynamics only, 



(22) 



Ts l + exp(-W^(t)/T) ' 

the total dissipation is zero, because in Eq. (|T^ the symmetric function of time (Tr{t) is 
multiplied by the antisymmetric function W{t). 

According to Eq. (|l^), the correction to the adiabatic limit, Aa,. = — <^r'^ satisfies the 
equation Act,,. = — (Xr/Fs; hence, as the first approximation in (3 we can use 



Aa, = -<7Fs. 

After the integration by parts of Eq. (|1^) the total dissipation is given by 

dW{t) 



(23) 



W{t) —r^ dt 



AcTr dt . 



(24) 



Combining Eqs. ( P^ ^ 
the total dissipation 



di J — oo dij 

(|T^), (|TH|), and (|^), we obtain the low-speed approximation for 



£ = n[3T = K 



ae 
GT 



(25) 



K = I ^ dx ~ 0.426 . 

-oo 



(26) 



x{l + e-^)3 

In this quasi-reversible regime £ <ti T; notice also that in this limit the total dissipation 
decreases when temperature increases. 

In the opposite limit, (3^1, the speed of energy change is so high that the tunneling 
occurs only at > and only once - see Fig. 0, without numerous back and forth processes. 
(Actually, this limit is the only one possible at T = 0.) Then the solution of Eq. ( [I6l) is 



ar{t) = 1 - exp (^- r+{t)dt^ 

( aGt'^\ 
= l-exp , t>0. 



(27) 
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and the average total dissipation 

S = (7r/5/2)^/2y = {'Re^al2Gf/^. (28) 

Hence, in the case (3^1 the average energy dissipation is much larger than (and independent 
of) the thermal energy, indicating the thermodynamic irreversibility of the switching process 
in this limit. In Fig. ^ the low-speed and high-speed approximations (Eqs. (^Sj) and (|28|)) 
are represented by dashed lines, while the solid line shows the result of the exact calculation 
using Eq. (0), which gives a natural crossover between these two limits. 

Equation (|^) is valid only for the case of the continuous spectrum of electron energies 
in the conducting islands. It is easy to consider a different case when the electron energy 
in the middle island is strongly quantized so that only one level is involved in tunneling, 
while the energy spectrum of the edge island is still continuous. (This situation is possible 
when the middle island is much smaller). Then the "orthodox" theory should be somewhat 
modified,&^^i however, for our purpose the only important change is that Eq. (^) should be 
replaced with 

TiW) = — 

T{W) + T{-W) = To = const, (29) 



where Fq is a constant which characterizes the tunnel barrier transparency. Equation (^0[) 
yields that the average total energy loss during one switching is given by a simple formula 

£ = a/To (30) 

for arbitrary switching speed (see Fig. |]). It is evident that quasi- reversibility {£ <^ T) is 
possible in this well, if the switching speed is low enough: a ^ T/Fq. 

In both cases considered above, electron energy relaxation in islands with continuous 
spectrum (implied by Eqs. (|^) and (^)) is the source of dissipation in the system. In the 
important case when electron energy is quantized in all islands this may not be true, so that 
this limit requires a completely different treatment. Our hope is to complete an analysis of 
this important case in the near future. 
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To conclude this section, notice that our model allows not only the total dissipation 8 
to be calculated, but also the time dynamics of the heat transfer between the system and 
the environment ("heat bath") during the switching process to be followed. Dashed lines 
in Figs. H and || show the results of such calculation using Eqs. (^ and (|29|), respectively. 
During the first part of the switching process (when W{t) < 0) the energy £i = —S{0) is (on 
the average) borrowed from the heat bath which, hence, is cooled. During the second part of 
the process (W(t) > 0) the average energy £2 = 8 — S{0) > £1 is returned back to the heat 
bath. Notice that in the adiabatic limit £1 = £2 = T \n2 independently of the exact model 
used for r(iy). This result follows from Eqs. (|19D and (|2^) after the integration by parts: 



Notice that this is valid for any gradual function W^(t), if only iy(0) = 0. 

The generality of this result is due to the direct relation £{t) = —TASmfit) between 
the energy and the informational ("Shannon") entropy Sinf of the degree of freedom used 
to code information (for a more detailed discussion, see Appendix B). In fact, in the instant 
when W = 0, the system may be in either of two states {am = (Xr = 1/2)) hence ASinf = ln2 
has been acquired in comparison with the definite initial state {am = 1,0-^ = 0). By the end 
of the switching the informational entropy is restored to the initial value since the state is 
definite again {am = 0, = 1). Finite switching speed decreases £1 and increases £2 (see 
the dotted lines in Figs. || and ^ leading to a positive total dissipation £ = £2 — £1. 



The general thermodynamic arguments lead to the conclusion that erasure of informa- 



Appendix B). During the switching OFF^ON of a cell in any SET-Parametron circuit, the 
amount of information is not changed, thus allowing arbitrary small energy dissipation in 
the slow-switching limit. However, for switching ON— s>OFF the lower bound on dissipation 
is determined by logical reversibility of a particular circuit. 




(31) 



VIII. REVERSIBLE COMPUTATION 



tion necessarily requires an energy dissipation of at least Tin 2 per bit (see Refs.E3E3 and 
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The SET Parametron shift register is obviously a logically reversible circuit, because 
during cell switching to OFF the information is preserved by the next cell. It is easy to 
check (see Fig. ^ that the sign of the energy difference A = ei — Er between two digital 
states does not change during the ON phase because of the back influence from the next 
cell. So the cell stays in the lower-energy state and the analysis of the energy dissipation 
during switching ON^OFF is equivalent to that of the switching OFF^ON (see Fig. |10|) . 
A similar small-dissipation case is realized in the SET Parametron fan-out circuit, because 
during the ON— i>OFF switching of the last cell of the input line (cell F in Fig. the proper 
sign of A is maintained by both first cells of the output lines. 

The situation is different for the NAND/NOR gate because any gate consisting of 
two inputs and one output is logically irreversible and, hence, has a lower bound for 
dissipation.&Ei For two uncorrelated input streams of bits with equal probabilities of unity 
and zero, the informational entropy before the logic operation is S^^j°^^ = — ln(l/4) while 
after the completion of the logic operation and erasure of the input information it becomes 
^inP^ = — ln(l/2). (Notice that the informational entropy decreases. This seems to con- 
tradict the apparent partial loss of information at computing. Actually the uncertainty is 
partially lost, not the information, because in the Shannon formalism we should treat the 
input data as unknown, and the data bit count is decreased by the gate.) The entropy 
difference determines the lower bound for the average dissipation S per logic operation 
S>T{Stj"'^-S!^}n=T\n2. 

Actually, in the SET Parametron realization shown in Fig. ^ the average energy dis- 
sipation is much larger than this lower bound. If two input bits are different, the energy 
difference A changes its sign during ON state of either cell A or B. In this case the energy 
dissipation during switching to OFF is comparable to |A| (see Fig. |10]) which should be 
much larger than the thermal energy because the condition |A| ^ T is necessary to ensure 
small error probability (see next Section). 

To realize reversible NOR and NAND operations (which would provide small dissipation) 
using SET Parametron cells, gates with two inputs and three outputs can be used (see Fig. 
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|TT]). This idea was sugg estedO for the Parametric Quantron logic gates and can be directly 
applied to the SET Parametron. The input information is copied to the first cells of two 
additional shift register lines. If their coupling to the last cells of input lines is stronger than 
input-output coupling, the proper sign of A is always maintained, and the energy dissipation 
is arbitrarily small in the slow switching limit. 



the SET Parametron cell switching (later we will briefly discuss also the "nonclassical" errors 
due to cotunnelingii). The error probability P is given by ai(oo). Let us first assume T = 
and calculate the "dynamic" error which occurs when the switching speed a is too high, 
and the system remains in the initial (symmetric) state up to the moment when tunneling 
to the upper energy level becomes possible (see Fig. |^ and Eqs. (0) - (|l^)). Since there 
is no tunneling back and forth at T = 0, the error probability can be found simply as the 
time integral of the rate of erroneous tunneling multiplied by the probability am that no 
tunneling has yet occurred: 



IX. BIT ERROR RATE 



Kinetic equation ( [T3| ) allows the calculation of the rate of "classical" digital errors during 




(32) 




(33) 



For the "orthodox" model of the tunneling rate given by Eq. (|^), 



Pdyn = i^7exp( 




(34) 
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(35) 
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where 7 = ae^ jGb?' . In order to keep Pdyn <^ 1, 7 should be much less than 1, then K 
Equation (Q) shows that the dynamic error decreases exponentially with the decrease of 
the switching speed and even faster than exponentially with the increase of A (factor in 
the exponent). 

For sufficiently small a and/or large A the main contribution to the error probability 
will be due to the thermally activated processes which populate the symmetrical state "m" 
during the passage of energy Sm across Si. The probability of this "thermal" error is given 
by the simple formula 

Ptherm = exp(-A/T) (36) 

and it prevails over Pdyn when T ^ ae^/GA. In the case when both errors are of the 
same order, the result can be found by the numerical solution of Eq. (|T3|). The total error 
probability can be estimated as the maximum of the two analytical results presented above: 

P ~ mayi{Pdyn, Ptherm)- 

If instead of the "orthodox" model we assume that only one discrete energy level of the 
middle island participates in tunneling, then using Eq. ( P^D for the tunneling rate we obtain 
from Eqs. ( |52D and (^) the following probability of the dynamic error (T = 0): 

Pdyn = ^ exp(-— Eo) . (37) 

The thermal error is still given by Eq. (j3^), and it prevails over the dynamic error if T 3> 
a/To. 



One more possible source of errors is the higher-order quantum process of cotunnelin; 




when two or more electrons tunnel simultaneously through different junctions. For illus- 
tration, the lowest energy diagram in Fig. |I]b shows the situation when the charge state 
with higher energy is occupied, and the digital information in the cell is preserved by the 
energy barrier (higher energy of the symmetric state) due to the bias field. According to 
the orthodox theory, single-electron tunneling in this case is impossible (at sufficiently low 
temperature). However, the second-order cotunneling, i.e. simultaneous tunneling of two 
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electrons (through both junction) brings the system into the lower energy state and, hence, 
is energetically allowed. This process changes the sign of the cell dipole moment and can 
lead to the digital error. 

This type of error can occur, for example, in the 3-phase shift register during the phase 
when the bit is stored by only one cell, and the long-range interaction with cells carrying 
other bits (nearest cells are OFF) provides uncontrolled sign of the energy difference between 
"/" and "r" states. The erroneous bit will then propagate along the shift register. 

Though a detailed analysis of cotunneling was not a goal of this work, we should notice 
that several means are readily available to reduce the resulting error probability. First, 
because the rate of m-electron cotunneling scales as {GRq)^ (where Rq ~ 6.5 kVt is the 
quantum unit of resistance) while the single-electron rate is proportional to the first power of 
G, the decrease of G decreases the relative importance of the cotunneling processes. Another, 
more powerful method is to increase the smallest order m of possible cotunneling processes. 
This can be done, for example, by increasing the number of islands per cell. If the cell 
consists of 5 islands and the ON state corresponds to an electron on one of the outer islands, 
then at least the 4-th order cotunneling (m = 4) is necessary to switch between logical zero 
and unity; a further increase of the number of islands per cell makes the cotunneling rate 
negligible even for relatively large tunnel conductance. Finally, it is possible to increase 
the minimal cotunneling order while still using three-island cells, by increasing the number 
of cells which store the same bit. This can be done, for example, by the increase of the 
number of phases in the operation of the shift register. In the realization shown in Fig. ^ 
this may easily be achieved by a decrease of the angle 9 between the planes of neighboring 
cells. If the bit is stored by k neighboring cells then the error can occur only if all these k 
cells will simultaneously change their polarizations and if the final state has a lower energy 
(any "partial" switching would cost at least one cell-cell interaction energy and, hence, is 
thermodynamically forbidden at least for not too large k). So, the lowest order of erroneous 
cotunneling is 2/c, and the linear increase in "hardware" allows the exponential reduction of 
the error probability. This method is applicable to any logic gates which use shift registers 
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as their input and output lines. 



X. DISCUSSION 

In this paper we have presented a functionally complete logic family based on the SET 
Parametron cells, including inverter, fan-out-two, NAND and NOR gates. All the circuits 
may operate correctly in a common range for the bias field amplitude. The important 
advantages of the SET Parametron logic are wireless operation and extremely low power 
dissipation possible in the quasi-reversible mode of operation, with the energy loss much less 
then ln2 per bit. 

While the realization of room temperature operation of the SET parametron logic re- 
quires sub-one-nanometer fabrication technology which is hardly available at the present 
time, the operation of simple SET Parametron circuits (with clock wires) can be readily 
demonstrated experimentally at lower temperatures. An obvious choice is the aluminum 
shadow-mask evaporation technology, widely accepted for singe-electron device fabrication 
(see, e.g., the coUectioni) . Figure |1^ shows the sketch of a possible layout of the two-cell SET 
Parametron shift register-inverter. (For the sake of simplicity, the artifact islands typical for 
this technology are not shown). Each cell consists of three islands. The capacitive gates are 
used for the application of rf "clock" fields which cause the switching processes, and simul- 
taneously for the application of dc fields to compensate random background charges. Input 
gates A and B (for the initial demonstration, one gate is sufficient) determine the polariza- 
tion of the first cell during its OFF-^ON switching. This polarization in turn determines 
the polarization of the next cell during its OFF— i>ON switching. The final polarization of 
the latter cell is sensed by a capacitively coupled single-electron transistor. (If a multilayer 
fabrication is available, this layout may be further improved to provide wider parameter 
margins: an overlap (without tunneling) of the edge islands of two cells would increase 
their capacitive coupling, while a similar overlap of the bias gates with the islands would 
allow smaller cross-talk.) The operation temperature of such a SET-Parametron circuit with 
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100-nm-scale tunnel junctions will be in sub-l-K range. 

Coming back to the wireless realization, let us estimate the parameters of a possible 
implementation of the device using the conducting (e.g., metallic) clusters as islands. For a 
cluster diameter 2R about 5 nm (which is at the limit of present-day direct e-beam writing 
techniques) the charging energy ~ e'^/SireeoR is about 0.15eV (where e ~ 2 is taken 
as a typical dielectric constant for the organic materials which are the natural candidates 
for the cluster environment). For the layouts considered in this paper the typical energy 
difference A between "/" and "r" states is about 0.2Ec leading to A ~ 0.03 eV. Requiring 
the probability of the thermal error to be less than 10^^'^ per switching, Eq. (|36D yields the 
maximal operation temperature of about 15 K. Assuming the same value for the dynamic 
error (then 7 ~ 0.025 - see Eq. (^4])) and taking into account the particular geometry (Fig. 

we obtain the maximum clock frequency i^max ~ 5x lO^^GA/e^ that corresponds to about 
10^ Hz for our parameters and G ~(10^ Ohm)~^ (higher G would make the cotunneling a 
problem). In this case the power dissipation (see Eq. (p8|)) is about 4(717/2) -"^/^z/A ~ z/A ~ 
5 X 10^^^ W per cell. This extremely low dissipation would make possible an integration 
level up to 10^^ cells per cm^ (i.e. 10^ nm^ per cell), limited by the cell size, since at 15 K the 
heat flux about IW/cm^ can be easily managed without the circuit overheating. To achieve 
the quasi-reversible mode of operation the frequency should be even lower: z/ ^ 5 x 10^ 
Hz for our set of parameters at T = 15 K. Accepting, for example, u = 10® Hz we get a 
power dissipation of only about 5 x 10^^^ W per cell. This figure indicates that as far as 
power dissipation is concerned, three-dimensional integration of SET Parametron circuits 
is quite feasible. Despite the not very spectacular clock frequency in this regime, the total 
computing performance of a 3D system can be very large. 

The largest problem with SET Parametron circuits (besides their fabrication) is that 
their operation requires well-defined background charges. The allowed fluctuations are only 
on the order of O.Ole (this number corresponds to the maximum tolerated angle of the 
cell tilt - see Section V). This is a common problem for any single-electron logic (see for 
example Refs.SSS). However, if /when a molecular technology becomes available for the 
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implementation of single-electron and other nanoscale devices, SET Parametron cells may 
be reproducible on the molecular level, and well defined background charges seem to be 
achievable in principle. 

From the estimates above we see that in order to increase the operation temperature of 
SET Parametron circuits up to 300 K or even up to 77 K, sub-nanometer-scale islands are 
necessary. In this case the energy level discreteness&il may result in radically new features 
of the SET Parametron operation, and should be taken into account at its quantitative 
analysis. Such analysis is outside the scope of the present paper. We would just like to 
notice in passing that the operation of SET Parametron circuits in that mode would be 
rather close to the recently proposed new versionil of the so-called "Quantum Cellular 
Automata" .0 This new version, put forward under the keywords of "adiabatic switching" 
and "pipelining", may remove the principal difficulties^ of the initial suggestion. We are 
not, however, aware of any detailed analysis of this new family of logic devices. 
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APPENDIX A: 

In this Appendix we discuss the method of multiple electrostatic images used to calculate 
the inverse capacitance matrix C^^ and vectors D (island dipole moments) which describe 
the infiuence of the external electric field for the arbitrary system of conducting spheres. 
This method was also used for numerical calculations in another work.§ 

It is well known that the electrostatic field of a point-like charge q located at a distance 
d from the center of an uncharged sphere having the radius r, may be treated as the net 
field of the charge q and a pair of fictitious charges located inside the sphere, in free space. 
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The image charge —qr/d is located at the distance r'^/d from the sphere center (in the 
same direction as the charge g), and the compensating charge qr/d is located at the sphere 
center. The total field of these three charges makes the sphere surface equipotential with 
= q/Aireod. 

Let us use this method to calculate the inverse capacitance matrix of two spheres with 
radii r and R, with distance d between the centers. For that, we need to calculate the 
electrostatic potential when a charge q is placed on the sphere R. First, let us create the 
pair ^gr/c? inside the sphere r (as if we had a point-like charge q). Three charges provide the 
equipotentiality of the sphere r, but not sphere R. So, we need to create two more charge 
pairs inside sphere R to restore its equipotentiality. One pair (image of the image charge) 
will have charges rkqd/r x R/{d — r'^/d) and the other pair (image of the compensating 
charge) will have charges ^qd/r x R/d (notice that one charge from each pair is positioned 
at the center of sphere R). Now the equipotentiality of the sphere r is lost again and we 
need to create four new pairs inside it (actually, three pairs because the position of two 
charges at the center of R coincides). We can continue this procedure until the charges of 
new pairs become sufficiently small, and the sphere potentials calculated at each iteration 
converge with the required accuracy. 

The same method is trivially generalized to an arbitrary number of spheres: we need to 
make images of all charges in all spheres. The method is very simple mathematically but 
cannot be applied in a straightforward way if the number N of spheres is large, because 
the number of required image pairs scales as A^" where n is the number of iterations. If 
N 100 and n ~ 20 (that is a typical number of necessary iterations if the spheres are close 
to each other, and we require a good accuracy), then the computer memory will obviously 
be insufficient. 

To solve this problem we use three ways (levels) of storage of information about the pairs, 
depending on their magnitude (the pair magnitude can be characterized by two numbers: 
the distance between the charges and the dipole moment of the pair). For "large" pairs 
we store the whole information (position and the charge magnitude). When the image of a 
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sufficiently "small" pair is calculated, we store only the location and the magnitude (dipole 
moment) of the dipole consisting of two image charges, and the total charge at the sphere 
center. The lowest level of the information representation is to consider the dipole being 
located at the center of the sphere, so we can sum up all the dipole moments. 

This modification of the algorithm allows us to calculate the capacitance matrix for a 
few hundred spheres with a typical distance between the spheres as small as one tenth of 
the sphere radius using a modest personal computer. (More closely located spheres require 
a larger number of image charges.) 

The calculation of the D vectors is analogous to the calculation of the inverse capacitance 
matrix. For example, to determine the influence of field we calculate ffist the dipole 
moments (which will have only x-component) produced by the field in the independent- 
sphere approximation. Then we use the iteration procedure described above to restore the 
equipotentiality of all spheres. At this stage image dipoles (with all spatial components) as 
well as image charge pairs appear. The island potentials for the unit external field are 
the components of the D vector corresponding to x-axis (we need three D vectors for three 
dimensions) . 



APPENDIX B: 

In this Appendix we discuss the heat transfer during the adiabatic transition (switching) 
between an arbitrary number of charge states. This analysis can be used for the fan-out 
circuits and for SET Parametron cells consisting of more than three islands. The formalism 
was developed long ago&B and constitutes the basis of the reversible computation analysis 
giving the lower bound for the energy dissipation for irreversible logical gates. 

The infinitesimal heat transfer to the thermal bath can be written as 

dQ = ^{Si - Ej) Ci Ti^j dt= - ^Ei dai , (Bl) 

where Ei is the (free) energy of the state i, (Ji is its probability, dat is the probability 
change during the interval dt, and Fj^j is the rate of transition from state i to state j (i.e. 
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the corresponding tunneling rate). In the thermodynamically reversible adiabatic limit the 
probabilities satisfy the equilibrium distribution 



(Ti = ai^ad = Z ^exp(-£i/T), Z = ^ exp(-£i/T). 

i 

Using the definition of Shannon's informational entropyH§ 
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and using the evident equation dci = we get 
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^(i(Ji(l + Ino-j) = T ^"^Eidai. 
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We see that in the adiabatic case 



dQ = -TdS, 



inf- 



(B5) 



It is not difficult to prove that in the general (nonadiabatic) case 5Q > —TdSinf- Let us 
introduce the difference 



Using the general master equation 
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to substitute rfcxj in Eq. ([B6|) , making the resulting expression symmetric over indices i 
and j, and using the general thermodynamic relation for the tunneling rates Ti^j/Tj^i = 
exp((£j — ej)/T) (c.f. Eqs. (plSf) and (^)) we obtain 
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After the simple transformations and using Eq. ( p2|) we finally get the expression 



[B8) 
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x{ajeM-h^r^)-(^i] (B9) 



which is obviously positive or zero. This proves the general inequality 



dQ>-TdSinf. (BIO) 

Notice that if we introduce also the "usual" entropy S so that SQ — T dS, then the total 
entropy St = Sinf + S is constant in the adiabatic case. In the general case, the total entropy 
is a non-decreasing function of time, dSt/dt > 0, while the decrease of S is not forbidden. 
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FIGURES 

FIG. 1. SET Parametron: (a) the basic cell and (b) its energy diagrams for three values of 
the bias field Ec (for discussion, see the text). 

FIG. 2. Phase diagrams of stable charge states of a SET Parametron cell with qio = 0, 
R/r = 1, d/r = 3, b/r = 1, where 2d is the distance between outer island centers while b is the 
y-axis shift of the middle island center - see Fig. |^a. (a) The cell is charged with one extra electron; 
(b) electrostatically neutral cell. AEc and AEs are the diagram periods for the clock and signal 
electric fields, respectively. Bistable regions are shaded. Thick line in (a) illustrates the switching 
from OFF to ON states, while rectangles in (b) illustrate two possible operation modes of the cell. 

FIG. 3. The general idea of the three-phase shift register using parametron-type cells. Each 
line shows the logical state of the shift register at a particular phase of the clock cycle, dashes 
indicating the monostable OFF state and I representing bistable ON states. 

FIG. 4. The top (left) and side (right) views of a shift register based on the single-exciton 
parametron cells operating in the symmetric mode. The clock field vector E{t) rotates in plane 
y — z. Digital bits are coded by the direction of the cell polarization and are propagated from the 
top of the figure to the bottom, over 6 cells during one clock period. 

FIG. 5. Time diagram illustrating the operation of the shift register shown in Fig. ^. Three sine 
curves show the in-plane component of the clock field for three neighboring cells. Dots correspond 
to the points where the middle curve enters or exits the bistable ON regions (shaded), so that the 
middle cell changes its charge state. 

FIG. 6. Circuit which can be used for the fan-out-two splitting of a signal (if the propagation 
direction is from the bottom to the top), or as a logical gate NAND or NOR (for opposite signal 
propagation direction). The asymmetry required for NAND and NOR gates can be created by 
the adjustment of the background charges on the edge islands of cell F, or by application of local 
electrostatic field (Es). 
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FIG. 7. Energy diagram of the OFF— >0N switching. The system, initiahy in state "m", is 
switched to state "r" (switching to "1" would give a digital error) . Numerous tunneling events (back 
and forth) occur in the slow switching case /? ^ 1 (a particular result of Monte-Carlo simulation 
for /? = 0.1 is shown), while for /? ^ 1 there is only one tunneling event. 

FIG. 8. Components of the energy exchange between the parametron and the heat bath as 
functions of the process speed a = dW/dt. Dotted lines: average energy flow £i from the heat 
bath to the parametron during the first part of the process (when W < 0), and the average flow £2 
from the device back into the heat bath during its second part (when W > 0). Solid line: the net 
energy dissipation £ = £2 — £\. Dashed lines show the low-speed (Eq. (p5|)) and high-speed (Eq. 
(p8|)) asymptotes. 

FIG. 9. The same as in Fig. ^, but for the case of the discrete energy spectrum of the middle 
island. 

FIG. 10. Energy diagram of the ON— >OFF switching. If the system initially occupies the 
state with lower energy, arbitrary small energy dissipation is possible at slow switching (simulation 
for /? = 0.1 is shown). On the contrary, if the switching starts from the upper state, the energy 
loss is finite and of the order of |A|. (For this case, Monte Carlo simulation has given two almost 
simultaneous tunneling events shown by arrows). 

FIG. 11. Schematics of logically reversible NAND/NOR gate with two input shift-register lines 
and three output lines. The input information is preserved in two additional output lines. 

FIG. 12. Sketch of a possible layout of the two-cell SET-Parametron shift register suitable for 
the standard double-angle shadow evaporation technique (the artifact islands are not shown). 
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